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ABSTRACT 

The lack of hydrogen in spectra of type la supernovae (SN la) is often seen as troublesome for single- 
degenerate (SD) progenitor models. We argue that, since continued accretion of angular momentum can prevent 
explosion of the white dwarf, it may be natural for the donor stars in SD progenitors of SN la to exhaust their 
envelopes and shrink rapidly before the explosion. This outcome seems most likely for SD SN la progenitors 
where mass-transfer begins from a giant donor star, and might extend to other SD systems. Not only is the 
amount of hydrogen left in such a system below the current detection limit, but the donor star is typically 
orders of magnitude smaller than its Roche lobe by the point when a SD SN la occurs, in which case attempts 
to observe collisions between SN shocks and giant donor stars seem unlikely to succeed. We consider the 
constraints on this model from the circumstellar structures seen in spectra of SN 2006X and suggest a novel 
explanation for the origin of this material. 

Subject headings: binaries: close — supernovae: general — white dwarfs 



1. INTRODUCTION 

The long debate over the progenitors of type la supernovae 
(SN la) has gained in urgency partly due to the use of SN la 
observations in cosmology, since it is somewhat embarrass- 
ing to use SN la for precision measurements when we don't 
properly understand how they are produced. The two broad 
classes of progenitor models are the 'single-degenerate' (SD) 
and 'double-degenerate' (DD) scenarios. In the former case, 
a carbon-oxygen white dwarf (CO WD) gains mass from a 
non-degenerate companio n until the WD reaches the point of 
explosive carbon burning (IWhelan & ibenl[T973l). In the lat- 
ter, the merger of two W Ds leads to the SN la dWebbinkl 19841 
llben & Tutuk"ov|[T984T) Pl 

Within each of those cateogories lies a variety of possibili- 
ties. The favoured SD systems include those with donors on 
the main sequence (MS) or the subgiant bra nch (collectively 
known as the supersoft c hannel ; Ivan d en Heuv el et al.l 

19941 

Lan ger et all 



1992; Rap paport et al 
iLan 



1997 



Li & van den Heuvel 
lHan & Po dsiadlowskil 120041 
Ivanova & Taaml 120041) and thos e with red-giant (RG ) 



2000; 



donors dHachisu etal 



19961 [19991 lHachisu & Katd [200lh . 



A third option, in wh ich the donor stars are helium-rich 
dlben & Tutukovl fl991l) has the advantage of avoiding hy- 
drogen contamin ation; however, it c annot produce a large 
fraction of SN la dWang & Hanll2010h . 

The known WDs closest to the Chandrasekhar limit are in 
SD systems. The WD in each of the novae U Sco, RS Oph, 
T CrB, and V445 Pup is inferred to have a mass very close 
to the Chandrase khar limitQ Furthermore, the observations of 
iPatat et al.l d2007l) suggest that the SN la 2006X came from a 
SD progenitor (see also section|4]i. There are also indications 
that the former donor stars in SD systems have been observed: 
in pop ulations of unusual WDs dHansenll2003l: lJustham et al J 
2009) and i n individual system s (Ruiz-Lapuente et al.l I2004t 
but see, e.g. JKerzendorf et al.ll2009l) . 

The strongest arguments against the SD channel are the 
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1 Clearly, diverse progenitors may help explain the diversity of SN la. 

2 Though we cannot tell whether the WD in any of these systems is a CO 
WD or an ONeMg WD. 



apparent lack of sufficient progenitor systems and the ob- 
servational lack of hydrogen in SN la spectra. Binary pop- 
ulation synthesis generally struggles to produce sufficient 
SD SN la to explain the empirical SN la ra t e (see, e.g., 
Han & Podsiadlowskil 120041 iMeng et alj [20091; iRuiter et alj 
2009, thou gh there are con siderable theoretical unc ertainties). 
bilfanov & Bogdanl d2010l) and lDi Stefanol d2010al) also argue 
that observed populations of SD SN la precursors cannot ac- 
count for all SN la, though this may be due to i ncorrect as- 
sumpti ons about the appearance of SD systems dDi Stefanol 
l20T0bh . 

The lack of hydrogen is a more acute problem, which could 
preclude almost any SD progenitors. Observations should be 
able to detect very small amounts of hy drogen (;$ 0.01M Q ) 
but none has been seen (Leonard 2007). Almost the en- 
tire envelope of a red-giant donor is b elieved to be strippe d 
by the shock wave from the explosion (Mar ietta et al.l 2000): 
such stripping should lead to significant contamination of the 
supernova spectrum by hydrogen. In the case of a main- 
sequence or subgiant donor, a smaller fraction of the enve- 
lope has been predicted to be stripped. Pak mor et al.l (2008) 
showed that the simulated amount of hydrogen stripping from 
main-sequence donors could still be within the observational 
limits, though only barely; this is u ncomfortable for the SD 
model (see also Meng & Yang 2010). Here we present a sce- 
nario in which the remaining envelope mass at explosion is 
naturally less t han the observational 0.0 IMq limit given by 
iLeonardl ((2007) and also very little of it can be stripped. As 
our model also significantly reduces the cross-section for in- 
teraction between the supernova shock and the donor star, it 
would unfortunately also reduce the opportunity for such in- 
teractions to be spotted in SN la light curves, as proposed by 
lKaser3 ( f2OT0l) . 

Section[2]explains how non-degenerate donors could finish 
mass transfer and shrink well inside their Roche lobes by the 
time of explosion, section [3] discusses whether it is plausible 
for this chain of events to occur in a significant fraction of SD 
SN la progenitors, and section |4] considers to what extent SN 
2006X constrains our model. 

2. EXPLODING AFTER ENVELOPE CONTRACTION 
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Table 1 

Changes in limiting WD mass with WD rotation 

100 



Situation 


Mass (M Q ) 


(a) Ignition of non-rotating CO WD 


1.378 


(b) Ignition of accreting WD with strong AM transport 


« 1.5 


(c) Max. for accreting WDs in Yoon & Langer model 


ss 2 


(d) Extreme stability limit of differentially rotating WD 


>4 



Note. — Masses are taken from: (a) Nomoto et al. ( 1984); (b) Uen- 
ishi, Nomoto & Hachisu (2003), Saio & Nomoto (2004); (c) Yoon & 
Langer (2005); (d) Ostriker & Bodenheimer (1968) 



Our mechanism for avoiding hydrogen contamination is 
built on two apparently robust effects. Firstly, accretion of 
angular momentum can increase the mass at which a CO WD 
is expected to explode and allow the mass transfer phase to 
continue for longer than is normally assumed. Secondly, if 
the donor star has a giant structure, then a natural end to the 
mass transfer phase occurs when the envelope mass becomes 
so small that the envelope rapidly shrinks. 

2.1. Postponing ignition during accretion of angular 
momentum 

Rotating WDs can avoid exploding or collapsing even if 
their masses significantly exceed the Chandrasekhar limit 
(Mch)- In the extreme case of differential rotation, a whit e 
dwarf is stable to 4 dOstriker & Bodenheimer! 119681) . 
More recently, lYoon & Langer! (120041 12005b found that SD 
SN la progenitors can reach rs 2 M Q before exploding. This 
may explain SN 2003fg, which appears to have h ad a super- 
Chandrasekhar WD progenitor dHowell et al. 2006)0 Even if 
accretion takes Mwd well above Mch, ignition may not occur 
until angular momentum is lost from or redistributed within 
the WD. 

In contra st to Yoon & Langer, ISaio & Nomotol (120041) and 

iPirol d2008) argue that angular momentum transport is very 
efficient, and hence that the accreting WD is unlikely to be 
far from solid-body rotation. How ever, solid body ro tation 
can also delay ignition : for example Uenishi et al.l (120031) and 
ISaio & Nomotol (120041) find that the ignition mass of the WD 
is raised by w 0. 1M Q even in the case of efficient angular mo- 
mentum transportFJ Table Q] summarises the WD mass limits 
for different models. We consider the effect of both the solid- 
body and differentially -rotating cases in section [3] 

There is an important caveat to invoking differential rota- 
tion and using the higher limit of w 2 M . Sustaining differ- 
ential rotation probably takes an accretion rate of > 10~ 7 M© 
(lYoon & Langerl l200l 12001 PI High enough mass trans- 
fer rates are far from guaranteed in all cases, e.g. for RG 
donors with low-mass cores. However, near-solid-body ro- 
tation seems sufficient for the scenario proposed here to be 
common & important (see section 3). 

2.2. Rapid contraction of giant envelopes 

If continued accretion can prevent the explosion, then the 
SN la will be postponed until after mass transfer ceases (and 
then until after sufficient angular momentum of the WD has 
been lost or redistributed). When the mass of a RG envelope 
(M env ) becomes too low for the star to maintain a red giant 

3 See also lChen & Uf2W% . 

4 In pure solid-body rotation, the WD stability limit is 1 AIMq . 

5 Of course, if differential rotation has been established, it continues until 
angular-momentum is redistributed. 
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Figure 1. The upper panel shows stellar radii as a function of envelope mass; 
the lower panel shows the thermal timescales of those envelopes. The labels 
next to the the curves show the fixed core masses used for each calculation 
and the curves in the lower panel become dotted when the stellar radius drops 
below 1 Rq. The knee in the radius-mass plot marks where the envelope 
starts to shrink; in all cases, this happens with M em below the O.OIM© limit 
for hydrogen detection given by Leonard (2007). 

structure then the envelope rapidly contracts on its thermal, or 
Kelvin-Helmholtz, timescale TKH,env, given by: 

GM star M env 

TKH,env « ^ (D 

where G is the gravitational constant, M stal and M env are the 
mass of the star and the envelope and R and L are the radius 
and luminosity of the star. We can scale that timescale to 
approximate stellar properties on the giant branch as: 

TKH,env _ / M star \ / M em \ t R \"' / L \"' 

48 yr ~ ^O.3M y ^10- 2 M Q / \1QR Q J \1Q 2 L J 

(2) 

Since the luminosity is a steep function of core mass, this 
timescale drops rapidly for donors higher up the giant branch. 
Figure Q] shows TKH,env when the envelope shrinks, calculated 
for a set of ^equilibrium stellar structures with different fixed 
core masses. We adopt a metallicity of 0.02 and use Eggle- 
ton's stellar evolution code (Eggletoa ll971tlPols et al.l ll995). 
In reality the envelope is depleted by shell burning from be- 
low as well as by mass loss from above; also the envelope 
will not be in strict equilibrium during mass transfer. For 
each core mass, then above a certain envelope mass the radius 
is roughly constant; below that envelope mass the envelope 
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shrinks rapidlyQ If the donor star in RS Oph is Roche-lobe 
filling , it has a core mass of « 0.4M Q (see, e.g. lJustham et al.l 
2009); based on these calculations, if the donor's envelope 
mass were to fall below ~2x 10~ 3 M Q , the envelope would 
contract in TKH.env < 1 yr. 

Once contraction of the envelope is complete then the ra- 
dius of the star is essentially insignificant; taking a separation 
of « 100^0 and a core size of rs O.IRq gives a fractional 
cross-section of ~ 10~ 6 for interaction between the supernova 
shock and the companion star. In addition, after contraction 
any remaining envelope will be tightly bound to the compact 
core. 

2.3. The delay before explosion 

After the end of accretion, the explosion will happen only 
after the WD redistributes or loses sufficient angular momen- 
tum (for differential or solid-body rotation, respectively). The 
timescale for this angular momentum loss or redistribution is 
uncertain. However, an upper limit to the angular momen- 
tum loss and redistribution timescales of 10 6 years has been 
claimed from mapping th e central density at ig nition to the 
expected nucleosynthesis ( Yoon & Langer 2005). But even if 
this delay is negligible then the finite time that carbon burning 
takes to runaway should be long enough to allow the envelope 
to contract. 

After the energy generation rate from carbon burning ex- 
ceeds the rate at which neutrinos can cool the core then con- 
vection currents are believed to postpone the explosion until 
the heating timescale becomes to o short compared to the con- 
vective turnover timescale (e.g. lArnettlir969l iLesaffre et al.l 
2006). The likely ~ 10 3 year delay due to this 'simmering' 
phase would allow time for exhausted RG envelopes to shrink, 
even without considering the time taken for redistribution or 
loss of angular momentum (see Fig. [TJ. 

2.4. Potential extension to the super-soft channel 

The above model naturally applies to SD progenitors of SN 
la with RG donor stars. Our proposed mechanism is signifi- 
cant even if it only applies to most such RG+WD SD systems. 
However, our scenario may also extend to at least some pro- 
genitors from the supersoft channel; such systems with late 
Case A or early Case B donors will evolve into systems with 
giant donors if there is enough mass remaining in the enve- 
lope, and if the WD does not explode first. If the envelope 
mass is already too small, then there is no problem to solve: 
the donor will shrink to become a WD or core-helium burning 
hot subdwarf in a wide binary (until, that is, the explosion of 
the WD breaks up the binary)Q 

2.5. One exceptional case 

A class of exceptions to our model might be systems con- 
taining donor stars which are able to reach helium ignition 
after growing their WDs to Mch whilst still posessing a sub- 
stantial hydrogen envelope. These donor stars would contract 
away from contact after igniting helium, allowing explosion 
to occur. Even in this case the donor star would no longer be 
Roche-lobe filling (by an order of magnitude or more). 

3. POPULATION EFFECTS 

6 Mass transfer is driven by the growth of the core due to nuclear burning. 

7 For the conditions which divide donors which will produce a WD from 
those which can ignite helium in their core see, e.g., Han et al. 1 2002). 



Despite the extended-accretion model for single degenerate 
progenitors described above, there exists a plausible explana- 
tion for why most SN la apparently explode within a narrow 
range of masses: the mass reservoir in the donor may usu- 
ally be too small to allow significantly super-Chandrasekhar 
masses. Below we suggest that most SN la might reasonably 
still explode with WD masses < 1.5M©. The WD mass dis- 
tribution at explosion directly relates to the question of what 
fraction of SD SN la are subject to our proposed mechanism. 

3.1. Principles & estimates 

The remaining donor mass at explosion depends on the 
donor mass at the start of the accretion phase (Mj ,,-), the ini- 
tial and final accretor masses (Mwd.i & Mwdj) and the mean 
accretion efficiency (i.e. the fraction of the mass lost by the 
donor that the WD manages to accrete and retain). There 
is considerable uncertainty in the population distribution of 
these: any model that could know them would have largely 
solved the question of which systems produce SN la. 

However, we can apply some constraints. The canonical 
limit for dynamically stable mass transfer limits the mass of 
a giant donor to be ;$ 1.2 times the mass of the accretor0 So 
for an initial COWD mass of < 1.0 Mq then for a RG donor 
any final core mass of 0.2 Mq or more would require a mean 
accretion efficiency above 50% before the final mass of the 
accretor could exceed « 1 .5 Mq (i.e. before we must appeal 
to strong differential rotation). 

Much lower - even negative - mass accretion efficiencies 
for such RG+ WD recurrent novae are common i n the liter- 
ature ( see, e.g., lYaron et all 120051) . Furthermore IChen & Lil 
(2009) find that rapidly rotating WDs can experience even 
lower accretion efficiencies than non-rotating WDs. 

For the supersoft channel, dynamical stability canonically 
restricts the mass of a radiative donor to be < 3 times the mass 
of the accretor, in which case overall accretion efficiencies ;$ 
20% would allow our model to be widely applicable in the 
restrictive case of solid-body rotations^ This seems challeng- 
ing without a s ignific ant and inefficent recurrent-nova phase. 
lHachisu et alj ([2008) have suggested that the limit for dy- 
namical stability could be significantly higher than expected, 
which would make it very difficult for the mechanism pre- 
sented here to apply to those SD systems. 

3.2. Population synthesis 

For a population of RG donor stars, tMeng & Yangj ([2010) 
calculated remaining envelope masses at explosion under the 
standard assumption that SN la occur at a CO WD mass of 
1.378 Mq. Extrapolating from their results, then for accre- 
tion efficiencies J$ 0.3, almost none of the systems in one of 
the two populations they synthesised would increase the WD 
mass by more than 0. 1M before the donor's envelope con- 
tracted. 

For more general SD systems. IChen & (12009) attempted 
to make super-Chandrasekhar SN la by considering accretion 
onto rotating WDs, but found 'in most cases the final masses 
of the white dwarfs are not significantly exceeding 1.4M Q '. 
They state that for an initial WD mass of 1.0M Q and for 
a metallicity of 0.02 then 'the maximum explosion mass of 

8 See, e.g., section 5 of lFfan et all J2002T) . 

9 For M(yo,-1.0 Mq, then Mjj S 3.0Mq; remaining a RG requires a 
final core mass for the donor of 5 0.5 Mq, leaving ;$ 2.5 Mq available 
to transfer, of which S 20% could be added to the WD whilst satisfying 
M WD,f g 1.5M . 
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white dwarfs is always less than 1.5 M \ Although many of 
their donors had substantial remaining envelopes, that work 
assumed that super-Chandrasekhar WDs exploded when the 
mass transfer rate dropped below 3 x 10~ 7 MQyr _1 , i.e. a crite- 
rion which estimates when differential rotation can no longer 
be sustained; future work should also consider that near-solid- 
body rotation can prevent explosion for WD masses 1 .5Mq . 

So it seems plausible that our mechanism could apply to 
most RG+WD SD systems, whilst an extension to the major- 
ity of other SD systems, though possible, appears more dif- 
ficult. There may be a balance between predicted SD rates 
and remaining envelope masses, since varying parameters to 
increase SD SN la rates (e.g. increasing mass accretion effi- 
ciencies) also seems likely to increase the incidence of donors 
with a significant remaining envelope mass at explosion. We 
encourage population synthesis studies to simultaneously in- 
vestigate the SD rates, distribution of SN la explosion masses 
and the fraction of hydrogen-contaminated SN la for a range 
of different assumptions. 

4. SN2006X 

The Na absorption features observed by iPatat et all (120071) 
in SN 2006X were interpreted as evidence of a SD progenitor. 
In particular, the distribution was reminiscent of a set of nova 
shells^] Si milar features have been seen in a small number of 
other SN la dBlondin et al.ll2009t iSimon et al.ll2009l 

4.1. Constraints 

The interpretation of a system with a RG donor star and 
recent nova outbursts constrains our model, as it provides an 
estimate of the time between the shrink i ng of the envelope 
and the explosion in this case. Patat et al] ( 120071) estimate that 
a time of ~ 25 years might have passed since the last nova 
outburst. This timescale for envelope contraction can easily 
be beaten (see Fig. [TJ, so at least there is no difficulty with 
proposing that the envelope has contracted since the last out- 
burst. However, expected simmering times are rather longer 
than this; it would be helpful for our model if the ^25 year 
timescale given by Patat et al. is an under-estimateOH Those 
few SN la that show these features may simply be fine-tuned 
or extreme cases. Perhaps, e.g., differen t ial rot ation trun- 
cates the simmering phase: Lesaffre et al. (2006) suggested 
that differential rotation could suppress the convective insta- 
bility of the core and hence potentially shorten the time before 
the thermonuclear runaway is triggered. 

4.2. A new model 

Alternatively, the interpretation as nova shells could be in- 
correct. A speculative but potentially elegant explanation oc- 
curs if, for some SN la - perhaps those where the donors 
manage to ascend highest up the giant branch - the remain- 
ing envelope is [partially] spontaneously ejected. For core 
masses <J O.4M , the dynamical and thermal times of the en- 
velopes become comparable as the envelope mass decreases 
(for the 0.4M Q core in Fig. Q] fjoienv ~ 0.2 yr at contraction, 
with f dyn w 0.04 yrH whilst for a 0.45 M core, we expect 
?KH,env < fdyiO- It has been argued that f K H,env < 10fd yn leads 

An alternative SD explanation was sugg ested by Hachisu et al. (2008). 

1 1 This timescale also constrains angular momentum redistribution or loss. 

12 For a total mass of O.4M , ! iyn « 0.08(i?/100R Q ) 3 / 2 yr. 

13 This inequality is true when our code fails before the envelope shrinks. 
Note, however, that the dominant causal timescale may be the sound-crossing 
time. 



to envelope ejectio n on the AG B, as it should lead to pul- 
sational instability (Soker 2008). In addition, RG envelopes 
approaching the tip of the giant branch can h ave binding ener - 
gies which are formally positive (see, e.g., lHan et alii 19941) . 
That combination of conditions makes an instability to pulsa- 
tional envelope ejection seem very plausible. The pulsational 
instability itself or shaping by the binary orbit might produce 
the discrete features observed by Patat et alPl 

5. SUMMARY AND CONCLUSIONS 

We have presented a scenario in which SD SN la could 
avoid having Roche-lobe-filling donors at the point of explo- 
sion. Rapid rotatior[3of the accreting WD stabilises it against 
ignition until after the H-rich envelope of the donor is almost 
exhausted and contracts. The companion at explosion is then 
a WD or hot subdwarf in a wide orbit. 

The key uncertainty is the fraction of SD SN la to which our 
proposed scenario applies. We have argued that it could plau- 
sibly work for the majority of systems in the RG+WD channel 
(for which the non-detection of hydrogen seems most prob- 
lematic), and perhaps for many other SD SN la (depending 
on, e.g., accretion efficiencies). 

This result s uggests that searches for hydrogen in SN la 
dLeona rd 2007) and attempts to see the co llision betwe en the 
SN shock and the donor stars in SN la (iKasenl 120101) may 
not distinguish between SD and DD progenitors. They may 
also be unab l e to d isfavour giant donors, as concluded by 
lHavden et alj (120101) . Hence it is especially important that 
we make the most of the information from those SN la which 
show evidence for circumstellar material. We have speculated 
on a novel origin for those features. 

We have also argued that, particularly if the absorption fea- 
tures seen in SN 2006X are due to nova shells, they could 
place constraints on the speed of angular momentum redistri- 
bution within, or loss from, the accreting WD and the duration 
of the simmering phase before explosion. Hence such obser- 
vations may not just help determine the progenitors of those 
SN la, but also given an insight into the physics operating in- 
side the doomed WD. 
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